Ultrastrong coupling between a cavity resonator 
and the cyclotron transition of a 2D electron gas in the case of integer filling factor 
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We investigate theoretically the coupling between a cavity resonator and the cyclotron transition 
of a two dimensional electron gas under an applied perpendicular magnetic field. We derive and 
diagonalize an effective quantum Hamiltonian describing the magnetopolariton excitations of the 
two dimensional electron gas for the case of integer filling factors. The limits of validity of the 
present approach are critically discussed. The dimensionless vacuum Rabi frequency ilo/'^o (i-e., 
normalized to the cyclotron frequency cjo) is shown to scale as yja nqw^, where a is the fine 
structure constant, nqw is the number of quantum wells and v is the filling factor in each well. 
We show that with realistic parameters of a high-mobility semiconductor two dimensional electron 
gas, the dimensionless coupling f2o/wo can be much larger than 1 in the case oi u ^ 1, the latter 
condition being typically realized for cyclotron transitions in the microwave range. Implications of 
such ultrastrong coupling regime are discussed. 
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The study of light-matter coupling in confined geome- 
tries has been in the last two decades a very interesting 
topic in atomic and condensed matter physics. The real- 
ization of extremely high quality mirrors and the ability 
to manipulate atomic states with extremely long lifetimes 
has made possible the observation of the strong light- 
matter coupling regime. In this regime, the coupling ex- 
ceeds the losses and it is therefore possible to spectro- 
scopically resolve the so-called vacuum Rabi splitting^'^. 
Strong light-matter coupling has been observed in vari- 
ous solid state systems, ranging from microcavity embed- 
ded quantum wellsM to superconducting circuits cou- 
pled to transmission line resonators^. In solid state sys- 
tems, thanks to the collective nature of the excitations, 
it is possible to go further the standard strong coupling 
regime and to enter a new regime, where the vacuum 
Rabi frequency is not only larger than the loss rate, but 
becomes comparable or larger than the bare frequency 
of the uncoupled excitations^"—. The ultrastrong cou- 
pling regime is interesting for example for non-adiabatic 
cavity QED phenomena reminiscent of the dynamical 
Casimir effectji^"— . Moreover, it can lead to a dra- 
matic modification of the quantum ground state (vac- 
uum) properties^ii^. 

In this paper, we investigate the coupling of the mag- 
netic cyclotron transition of a two-dimensional electron 
gas (2DEG) to the quantum field of a cavity resonator. 
We show that the dimensionless vacuum Rabi frequency 
^.n/uifj can be largely enhanced with respect to the case 
of intersubband transitions in a 2DEG without magnetic 
field. In particular, this is the case in the regime of 
high filling factors, obtained with relatively weak mag- 
netic fields and large electron densities, which can be 
obtained in state-of-the-art high mobility 2DEGs. We 
derive the second quantized Hamiltonian for such a sys- 
tem in the case of integer filling factor v and derive an 
effective Hamiltonian describing the cavity magnetopo- 
lariton excitations of the 2DEG. We show that since 
rifl/wo 3> 1, the diamagnetic A^-term of the quantum 



light-matter coupling Hamiltonian becomes dominant in 
such a system. The present work is relevant not only 
for the fundamental quantum electrodynamical proper- 
ties of a 2DEG in a unconventional regime. It may have 
important implications also in the low-frequency magne- 
totransport properties of the 2DEG embedded in a cavity 
resonator. 

The paper is structured as follows. In Sec. HI we in- 
troduce in detail the system and show why the mag- 
netic cyclotron transition can be ultrastrongly coupled 
to the vacuum field of a cavity. In Sec. [Hi we present 
the second quantized quantum light-matter Hamiltonian 
for the system (details about the derivation are given in 
the Appendix). In Sec. IIIIl we diagonalize an effective 
bosonic Hamiltonian to describe the magnetopolariton 
excitations of the 2DEG and show the resulting mode 
dispersions. Conclusions and perspectives are drawn in 
Sec. IVl 



I. PHYSICAL SYSTEM AND SCALING OF 
COUPLING 

We will consider a system consisting of multiple doped 
semiconductor quantum wells (QWs) in presence of a 
magnetic field B along the z axis (perpendicular to the 
QW plane). The QWs are embedded in a wire- like cav- 
ity resonator, as depicted in Fig. [U that confines the 
electromagnetic modes along two directions (z and y). 

In presence of a magnetic field, the electrons occupy 
highly degenerate bands (see Fig. [5]), the well-known 
Landau levels (LLs), separated by the cyclotron en- 
ergy equal to hwQ, where ujq = eB/{m*c) is the cy- 
clotron frequency (m* « 0.068 toq is the effective electron 
mass of the conduction band in GaAs). The magnetic 
length, associated to the fundamental Landau level, is 
Iq — \J h/{m*u]Q), while the degeneracy of each LL (tak- 
ing into account for the electron spin) is N = S/[-kIq), 
where S is the surface of the sample. 



Here, we will consider the case of an integer filling fac- 
tor v, i.e., electrons fill completely the first ly LLs (the 
LL index goes from n = Oton = i/— 1). Consequently, 
1^ = P2DEGS/Af where P2DEG is the electron density per 
unit of area in each QW. In this section, we wish to show 
a back-of-thc-envelope calculation of the vacuum Rabi 
frequency for such a system and compare it to the bare 
cyclotron transition frequency wq. This simple calcula- 
tion will allow the reader to grasp the essential quan- 
titative features and key parameters. A complete and 
rigorous derivation will be presented in Sec. and in 
the Appendix. As we are interested in the observation of 
high filling factor LLs, we have to consider a sample at 
cryogenic temperature. In the rest of the manuscript we 
will thus perform calculations at T = 0, even if our results 
are expected to be robust while the thermal energy does 
not exceed the cyclotron transition energy {ksT < hwo). 

Due to Pauli blocking and harmonic oscillator selec- 
tion rules, only electrons populating the level n — v ~ 1 
will participate to the light-matter coupling. Given that 
the typical cyclotron radius for an electron in LL n is 
r ^ loy/n, the electric dipole associated to the tran- 
sition from level n = v — 1 to level n = v will be 
d ~ e^ov^- Hence, larger filling factors produce a larger 
dipole (in a way reminiscent to Rydberg atomic states 
with large orbital principal quantum number). However, 
if N2DEG — P2DEgS is the number of electrons in the 
2DEG, only the number N^degI^ in the level n = — 1 
is optically active. It is known that in the presence of 
a collection of identical dipoles, the collective vacuum 
Rabi frequency is proportional to the square root of the 
number of dipoles^^. Hence, the collective excitation of 
the 2DEG enhances the light-matter coupling by a factor 

N2deg/v. If riQw quantum wells are identically cou- 
pled to the field, an additional enhancement of yjnqw 
can be obtained (we consider here that the QWs are not 
electronically coupled; this is therefore equivalent to sim- 
ply increase the density of the 2DEG). 

The vacuum Rabi frequency is therefore Mires = d x 
Eyac X \/ N2D EG n-Qw I v., where E^ac is the electric field 
associated to the cavity vacuum fluctuations. We will 
call V — L^S the volume of the cavity, with being 
the direction orthonormal to the 2DEG and S = L^Ly 
the surface of the 2DEG plane (cf figure [T]). Now, if we 
consider a cavity mode resonant with the cyclotron tran- 
sition, we have Eyac ^ \J ^77%; where e is the dielectric 

constant of the cavity spacer ( e ~ 13 for GaAs struc- 
tures). Finally, we find 




If we consider a half- wavelength cavity, i.e., Lz = Ao/2 
and ojo = 1^^' given that P2DEG'^^) = we obtain 

~ y/a nqwv, (2) 
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FIG. 1: Sketch of a cavity resonator embedding uqw identi- 
cal quantum wells (QWs), each containing a two-dimensional 
electron gas (2DEG) tn the xy plane. An uniform and static 
magnetic field B is applied along the z axis. 
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where a = |- ~ 7^ is the fine structure constant. Hence, 
he 137 ' 

the dimensionless vacuum Rabi coupling depends on the 
small fundamental constant a, on the number of QWs 
and on the filling factor in each QW. For very large fill- 
ing factors and a large number of QWs it is thus possi- 
ble to have ^ 1. For a given density of electrons, 

since ly cx we can thus always increase the cou- 

pling by lowering the magnetic field intensity, since for 
B ^ 0, V ^ -\-oo. One the other hand, for a given 
magnetic field, the filling factor increases with increasing 
the electron density, and so does the coupling strength. 
Of course, this description makes sense only if the cy- 
clotron resonance is well resolved (i.e., not quenched by 
the broadening). For example, as shown by the experi- 
mental work in Ref. 19, in GaAs 2DEG with a relatively 
high mobility {p = 1.6 10^ cm^V^^s^^), it is possible to 
have at the same time very high filling factors {ly > 100) 
with a well resolved cyclotron resonance for a magnetic 
field B = 0.01 T and a transition frequency in the mi- 
crowave range {fo — ^ ~ 30 GHz). 



II. QUANTUM HAMILTONIAN 

A. General considerations 

In condensed matter systems consisting of a collection 
of two-level systems, in the low excitation limit, it is gen- 
erally possible to describe the collective excitations as 
bosons. Starting from the Hamiltonian in terms of the 
electronic fermion operators and of the photonic boson 
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FIG. 2: Sketch of Landau levels with an integer filling factor v 
(ly ^ 1 is the regime considered in this paper). The cyclotron 
transition between levels v — 1 and v is resonantly coupled to 
cavity mode quantum field. 



operators (see Appendix), we have derived the effective 
bosonic Hamihonian describing the magnetopolariton ex- 
citations. In this case, the Hamihonian is exactly solvable 
by mean of an Hopfield-Bogoliubov diagonalization. This 
procedure has been already applied to other systems: 
for example, exciton polaritons in bulk materials^S and 
quantum wells^ii^, intersubband polaritons^- and polari- 
ton in atomic gases2^. Among the different possible ge- 
ometries, we chose to treat here the case of a photonic 
asymmetric wire {L^ 3> 3> Ly). This will permit us 
to treat the continuous dispersion along the x axis, while 
keeping only few discrete modes along the z axis and ne- 
glecting completely all the modes along the y axis except 
the lowest one, whose intensity profile is constant over the 
sample surface. While this choice is motivated both by 
the simplicity of the resulting Hamihonian and by con- 
siderations of experimental feasibility (using transmission 
line resonators) the extension to other one dimensional 
geometries does not present any particular difficulty. For 
sake of simplicity, we will consider the quantum wells to 
be placed in the middle of the photonic wire, at z = Lz/2. 

The reader familiar with the literature on polaritons 
will notice that here, even in the resonant case (one pho- 
tonic mode resonant with the cyclotron frequency) we 
will not neglect the higher lying photonic modes, as often 
done in the literature. We will instead formally consider 
all of them, cutting then the resulting infinite matrix 
in order to retain enough modes to be at convergence. 
While considering only one photonic mode remarkably 
simplifies the algebra, it is not a valid approximation in 
the system under consideration. The spacing between 
photonic branches being constant, in the resonant case, 
the condition flo/ojo > 1 implies that also the higher 
photonic modes are coupled to the transition. 



B. Coulomb interactions 

As well known in the context of the fractional quan- 
tum Hall effect, the Coulomb interaction can play a 
crucial role in two-dimensional systems of electrons un- 
der magnetic field. It is now well understood that the 
role of interactions is completely different depending on 
whether one considers the case of integer or fractional 
filling factor. In the case of fractional filling, correla- 
tion can become crucial, because a rearrangement of a 
many-electron configuration within a partially filled Lan- 
dau band does not cost any kinetic energy. As it has been 
said in the previous sections, we consider here only the 
case of an integer filling factor v in the peculiar weak 
magnetic field regime such that ly ^ 1. In this weak 
magnetic field limit, we have e^/(e/o) 3> huio: hence, one 
might a priori expect a significant impact of Coulomb 
interactions, which could produce for example a mix- 
ing of the Landau levels. However, it has been proved^i 
that screening leads to a renormalized Coulomb poten- 
tial which abruptly drops at the distance around two cy- 
clotron radii. Moreover, the effective interaction is much 
smaller than Hujq and this allows us to treat only the 
electrons belonging to level n = v — 1 (which is the 
only one optically active). In addition, Kohn's theorem- 
states that the cyclotron resonance is not affected at all 
by electron-electron interactions as far as we consider a 
translationally invariant system. This holds in the ab- 
sence of disorder and for a photonic wavevector q = 0. 
In the geometry we consider (cf sec. Ill Al and III Cp . the 
photonic wave vector always satisfies the condition that 
qlo <C 1 (even far from resonance). Indeed, since for 

3> 1 the vacuum Rabi frequency Qq can be even much 
larger than the cyclotron frequency ujq, the light-matter 
interaction in the cavity system appears to be by far the 
most dominant interaction. 



C. Cavity quantized electromagnetic field 

In the considered geometry, the vector potential can 
be written as 



A{r) = Aoir) + Ae™(r), 



(3) 



where Ao{r) is the applied uniform magnetic field di- 
rected along the z-direction (Aq — —ByUx in the Lan- 
dau gauge that we will use in this manuscript) and Aem is 
the contribution associated to the cavity quantum field. 
We introduce the photon wavevector : 



(4) 



where Qx can vary continuously, while qy and q^ are quan- 
tized {riy and are integer values). 

In the following, we will take the length such that 
the mode corresponding to q^ = 0^ qy = and qz = 
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is close to resonance with the cyclotron transition. As 
explained in Section [II Al the condition Ly « allows 
us to neglect all the modes with qy ^ Q. To simplify 
the notation, we will thus omit the qy = index. The 
electromagnetic vector potential is then written as : 




.nM. 



„T JJ* 



(5) 

where the operator Uq^^n^ is the bosonic annihilation 
operator for a photon belonging to the mode labeled by 
the wavevector {qx,qz = t^}- The spatial shape tlq^,n^ 
of the modes (the derivation is presented in Appendix E)) 
is given by : 



U 



( 



\ 



sin(^z) 
V ) 



(6) 



Being our quantum wells positioned in the middle of 
the photonic cavity, at z — Lz/2, only the odd photonic 
modes will be coupled to the electron gas. Moreover we 
introduce the matrix notation : 



a„ 



(7) 



ttq^ is thus a vector containing the photon annihilation 
operators for all the modes with different values of = 



degeneracy. We also neglect the Zeeman splitting of the 
LLs (the cyclotron transition frequency is the same for 
both spin channels). 

Analogously to the case of intersubband transitions^^, 
the bright mode creation operator associated to the cy- 
clotron transition is given by 



" nQwP2DEGJ 



So) 



(8) 



The prefactor is a normalization constant chosen such 
that, in the dilute regime under consideration, these oper- 
ators are approximately bosonic, [&g^,&J/] — Sq^,q'^- The 
bright mode is the collective excitation of the two dimen- 
sional electron gas which is directly coupled to the cavity 
photon mode. 

After some calculations that are detailed in Appendix 
[Cl we get the following Hamiltonian 



H -H^Landau ^ Hint ^" -^dia ^" -^cavityt (9) 



where 



H Landau — ^ ^^^0 ^(jx ) 



(10) 



Hint — 



D. Light-matter coupled system 

In the Landau gauge each electron state is indexed by 
three quantum numbers n, k and m (see Appendix [B] for 
details), n indexes the different Landau levels, k is the 
momentum component along x and labels the different 
electrons in each Landau level, m indexes the different 
subbands in the quantum well. Actually we show in Ap- 
pendix [C] that given the condition Lqw ^ L^, we have 
the selection rule m = m' , meaning that we can safely 
neglect intersubband transitions and thus drop the m in- 
dex. Beside, one could argue that all the energies in play 
are much smaller than the intersubband gap which rein- 
forces this approximation. 

In order to write the second quantized Hamiltonian 
for the coupled light-matter system, we will thus intro- 
duce the fermionic operator c|j^)J, that creates an elec- 
tron with quantum number k in the n*^ Landau level in 
the j*'' quantum well. We omit the spin quantum num- 
ber because the optical transition conserves the electron 
spin, hence both spin channels contribute equally to the 
light-matter coupling: it will simply appears as an added 



Hdia = "^a^^Dq^a^q^+a^^Dq^al^ 

<?x 

-|- al^ Da tta + a'J' Da a} , 

' Qx yx yx Qx yx — <?x' 

and 

Hcavity = ^ al^fujq^aq^, (11) 

Sx 

which represents the energy of the cavity quantum elec- 
tromagnetic field. Note that we have omitted the zero- 
point energy that does not play any role here and we have 
defined 

Uq^ = diag [ujq^,l U!q^^2 ' ' ' ' ' ' ] (12) 

as the diagonal matrix containing the photonic mode en- 
ergies. H^andau IS the energy of the collective excita- 
tion electronic excitations. The vector flq^ and the ma- 
trix Dq^ contain the coupling constants corresponding to 
Hint and Hdia respectively. J7g^ is the collective vacuum 
Rabi frequency and comes from the resonant coupling 
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between the cavity photons and the 2DEG. Hdia (dia- 
magnetic term by analogy with atomic physics) comes 
from the squared vector potential A^^. In Appendix 
[Cl we present the detailed derivation of these coupling 
constants. The final result is 



1 2TTe'^uJo nQwP2DEG 



e m* L, 



(13) 



D 



r2„ nj, 

q^ 
LOO 



with 



diag 



^ n 



9x,7' ■ 



(14) 

and the zeros at the even positions in Eq. ((T?| are 
given by the fact that, as apparent in Eq. ([6]), only odd 
photonic modes are coupled to the electron gas. 

If we consider only the resonant photonic mode, with 
frequency ujq^^n^ = wq, then the corresponding normal- 
ized vacuum Rabi frequency reads 




FIG. 3: The dimensionless vacuum Rabi coupling Slres/i^o 
versus the filling factor v for different values of the number 
nqw of quantum wells. We consider only the photonic mode 
{qx = 0, riz = 1} at resonance with the transition. Precisely, 
the resonance is defined such that u)q^^o,n^^\ = cjq. Other 
parameters are given in the text. 



1 2a riQwi^ 



(15) 



Apart from a geometric form factor of the order of 1, we 



see that 



— ^ oc ,/c 

Ll>0 V 



V as in Eq. ([2]) , obtained with a 



back-of-the-envelope scaling calculation. 



From Eg. (113)) it is clear that, for 



^ 1, we have 



Ores ^ ^res- lu the ultrastrong coupling regime, we 
conclude that the term becomes thus dominant over 
the vacuum Rabi coupling term. 

Because this Hamiltonian is quadratic in terms of a^^ 
and hq^ operators, it can be exactly diagonalized using a 
generalized Hopfield transformation^^. We introduce the 
normal modes ]}ql and •p'qf^ (magnetopolaritons) defined 
as 



where ftwg^ correspond to the magnetopolariton ener- 
gies, v{l is the vector {w^} , X^q^ , Y^f^ , and Mq^ 
is the infinite Hopfield matrix that is given by 



/ 



Ma 



Wo 

'9x *^<?x 



2D„ 



-2D„ iVl„ 



Jq^ - 2Dq^ iftq^ 



\ 







-Wo / 



(18) 



P. 



^qx "-gx 



Given the bosonicity of a^^ and bq^ operators, also the 
magnetopolariton operators satisfy the Bose commuta- 



tion rule 



(i) Hi) 
PqJ^Pql 



Si.i'S, 



9x,gi 



where the index i runs 



over all the polariton branches. The condition for the 
total Hamiltonian to be diagonal in terms of magnetopo- 
lariton operators is 



and the eigenvalue problem takes the form 



(16) 



(17) 



III. RESULTS 

By diagonalizing the matrix in Eq. (jl8p . we are now 
able to calculate the magnetopolariton energy disper- 
sions in the multimode case. For sake of simplicity, we 
have taken the lower photonic mode to be resonant with 
the cyclotron transition. In Fig. [31 we present results 
of the dimensionless vacuum Rabi frequency VLresl'^o as 
a function of the filling factor v for different values of 
the number nqw of quantum wells. It is important to 
point out that the dimensionless coupling depends only 
on riQwv- Hence, different values of the two dimen- 
sional electron gas density P2DEG, magnetic field B, and 
semiconductor effective mass m*, can give rise to the 
same normalized vacuum Rabi coupling flj-es/^O: pro- 
vided that the filling factor v stays constant. For exam- 
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pie, a filling factor oi = 100 with a number of QWs 
of ngvK = 100 gives a dimensionless vacuum Rabi cou- 
pling that takes the impressive value V.res/^Q — 3.6. This 
kind of parameters can be obtained using GaAs high- 
mobility samples with cyclotron frequencies in the mi- 
crowave range, as already discussed at the end of Sec. 

m 
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FIG. 4: (Color online). Thick solid colored lines: normal- 
ized frequency of the first four magnetopolanton branches as 
a function of the 2D electron gas density P2deg in the case 
where the energy of the first cavity mode {qx = 0, = 1} 
is resonant to the cyclotron transition. Thin black solid lines 
depict the energies of the photonic branches, which are not 
coupled to the cyclotron transition (for the selection rules, 
see description in the text, in particular Eqs. U3}) and 
Iil4\ l)- The two vertical straight lines indicate the densities 
p2DEG = 2 lO^^cm"^ and P2deg = 2 lO^^cm"^ (giving 
a dimensionless vacuum Rabi coupling Hres/i^o — 0.8 and 
Qres/i-^^o — 2.5 respectively): the corresponding frequency dis- 
persions as a function of the wavevector q^ are plotted in Figs. 
[3| andO Parameters : m* = 0.068mo (GaAs effective mass); 
B = 40mT (giving a cyclotron frequency ojq = lOOGHz/rad^; 
ngw ~ 50 and Lz ~ 0.25 cm which is chosen in such a way 
that 0)5^=0, nj=i ~ 1^0. With these parameters, a good numer- 
ical convergence is obtained for a cutoff number of photonic 
modes A > 7. 




FIG. 5: (Color online). Thick solid colored lines: normal- 
ized frequencies of the first four normal mode branches as a 
function of the wavevector q^ for P2deg = 2 lO^^cm"^ and 
B = 40mT which correspond to a filling factor v ~ 10. The 
dashed color lines represent the dispersions of the photonic 
modes with odd values of that are coupled to the cyclotron 
transition. The red dashed-dotted line depicts the bare cy- 
clotron frequency; the thin black solid lines represent the pho- 
tonic modes with even values ofn^, which are not coupled to 
the electronic system. Other parameters are the same as in 
Fig.\J\ 




In Figure m the thick solid lines represent the normal- 
ized energies of the first four magnetopolariton branches 
as a function of the two dimensional electron gas density 
P2DEG in the case where the energy of the first cavity 
mode {qx = 0, — 1} is equal to the cyclotron transition 
energy. Note that the thin solid lines depict the energies 
of the photonic modes that, due to the considered geom- 
etry, are not coupled to the cyclotron transition. Fig. [5] 
and [6] depict a typical dispersion of the normalized fre- 
quencies of the first four polariton branches as a function 
of the wavevector for densities P2DEG = 2 lO^^cm^^ 
and P2DEG = 2 lO^^cm"^ respectively. 



0^ < < 1 < < 1 < 1 

5 10 15 20 25 30 35 40 

q^(cm"'') 

FIG. 6: Same as Fig. [5l but with an electron density p2DEG = 
2 lO^^cm"^ (corresponding to a filling factor v ~ 100^. 



It is worth pointing out that for very large couplings 
the frequency of the lower branch asymptotically tends to 
0. This suggests that this kind of cavity excitations could 
affect significantly the magnetotransport properties of a 
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cavity embedded two dimensional electron gas in the low 
frequency regime. 



IV. CONCLUSIONS AND PERSPECTIVES 

In conclusion, we have presented a quantum model de- 
scribing the ultrastrong coupling between a cavity res- 
onator and the cyclotron transition of a two dimensional 
electron gas. The present approach holds for the case of 
integer filling factors. We have determined the second 
quantization light-matter Hamiltonian in terms of the 
electronic fermionic operators and of the cavity photon 
operators. We have derived and diagonalized an effective 
bosonic Hopfield-like Hamiltonian describing the magne- 
topolariton excitations. The ultrastrong coupling regime 
characterized by a vacuum Rabi frequency much larger 
than the cyclotron transition frequency can be achieved 
with high filling factors, which are compatible with state- 
of-the-art GaAs high mobility two dimensional electron 
gas and cyclotron transitions in the microwave range. In 
our present approach, we have not considered the case 
of non-integer filling factors: the nature of the magne- 
topolariton excitations may be qualitatively and quan- 
titatively affected by a partially filled Landau band, an 
issue that needs to be explored in the future. 

The role of Coulomb interaction between carriers and 
basic approximations have been critically discussed in 
Sec. (HTBI) . 

In the future, it will be interesting to explore the im- 
pact of the ultrastrong coupling on the magnetotrans- 
port of properties of a cavity embedded two dimensional 
electron gas. In fact, the energy of the lower magne- 
topolariton branch can become much smaller than hujQ 
for r^res/i^o ^ 1: meaning that even the low-frequency 
transport can be affected by the coupling to the cavity 
vacuum field. In particular, the Shubnikov de Haas os- 
cillations, the nonlinear response to applied microwave 
fleld^^"— , could be strongly influenced by the presence of 
the cavity. 

We would like to thank S. Barbieri, Y. Gallais, C. Sir- 
tori and Y. Todorov for discussions. 



to the mode labeled by the wavevector : 



sm f cos ( 
= q I sin 9 sin <; 
cos 9 




(A2) 



where Qx can vary continuously, while Qy and Qz are quan- 
tized (uy and Uz are integer values). The spatial shape 
Uq jj of the modes is given by2S 



/ ism{qzz)sm{qyy) 



sin((j2z) cos{qyy) cos 9 sin 
\ — cos{qzz) siii{qyy) sin9 j 



(A3) 



/ -isin(q2z)sin(gj^2/) sin(/)\ 
Wq,2 = Ne"'"^ sm{qzz) coa{qyy) cos(f> 

V / 

The normalization constant N reads 



N 



if qy = or qz 







(A4) 



otherwise. 



As it has been written in Sec. Ill Al and III CI the con- 
dition Ly << Lz allows us to neglect all the modes with 
qy 0. From Eq. (jA2l) . we see that qy = implies = 
and consequently only the polarization 77 = 2 is present. 
Omiting the qy — and ry = 2 indexes, we can see from 
Eqs. (Ml . (TMj) and (K3\i that the electromagnetic vec- 
tor potential and the spatial shape of the cavity modes 
are given by Eq. ([5]) and (jl]) respectively. 



Appendix B: Landau Levels in the Landau gauge 

The one electron wavefunctions in the Landau gauge 
{Aq = -Byux), reads 



Appendix A: Spatial shape of the cavity modes 



VnM.m{r) = Xk{x) (t>„iy - yo) Cmiz), (Bl) 



The general form for the electromagnetic vector poten- 
tial is 



where 



where the operator ag ^j is the bosonic annihilation oper- 
ator for a photon with polarization rj € {1,2} belonging 



Xk{x) 



i>niy-yo) = 



Akx 



II — 7/„ \ (y-yp) 



2/0 = klg is the so-called guiding center position de- 
pending on fc, ^m{z) is the confinement wavefunction of 
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the m-th conduction subband of the quantum well and 
Hn is the Hermite polynomial of degree n. With the 
chosen gauge the wavefunction is thus factorized along 
the three axis in a plane wave, an harmonic oscillator 
wavefunction and a confinement-dependent function. For 
an infinitely deep quantum well with width Lqw, we 
can find an analytic form for the confinement-dependent 
function 



writing A{r) = A^lr) + Aem{r) (as in Eq. ([3])), we can 
identify four different terms in the Hamiltonian, namely 



H 



Hi, 



H, 



dia 



H, 



cavity ■ 



(C2) 



■ cos 



odd 



(B2) 



As explained in section III D| we omit both the spin 
quantum number, the Zeeman splitting of the LLs and 
the m index that do not play any role here. 

Knowing that k = iimx/ with Ux € N, the many 
body electronic ground state reads 



i^)-nnn 



„0")t 



|0), 



(B3) 



-1 n— — 1 



H Landau — ^hlndau ^^'^ Hint — l^j ^inf H cavity 

describes the free quantum cavity electromagnetic field 
and has been expressed in terms of second quantized 
photon operators in Eq. (ITT|) . Hdia comes from the 



tU) 



term. H 



Landau 



is the Hamiltonian describing the elec- 



trons in the j*'' quantum well, interacting with the static 
magnetic field, giving rise to the Landau levels. Finally 
^int describes the interaction of the electrons in the j*'' 
quantum well with the cavity electromagnetic field. In 



the following we will express H^ll^dau^ ^ini and Hdia in 
terms of the a^^ and cl^\ operators introduced in section 
HTCl HTdI and Appendix E] 

^Landau Construction a diagonal operator, be- 

cause its eigenmodes are chosen as a basis for the transi- 
tion in second quantization formalism, namely: 



rU) 



where |0) is the empty conduction band state and c^^^^ 
is the fermionic operator creating an electron in the n*'' 
Landau level with wavevector k along x in the j*'' quan- 
tum well. 



(C3) 



Appendix C: Second quantization hamiltonian 

The minimal coupling Hamiltonian describing the elec- 
trons coupled to the electromagnetic field reads: 



cavity ; 



where the sum runs over all the electrons (in all the 
quantum wells) and Hcavity describes the free electro- 
magnetic field. Developing the square in Eq. (jCl|) and 



As already mentioned in Sec. IH Dl because the elec- 
tronic excitations are collective, we have to express 



H 



U) 



idau 



in the subspace of the collective excitations. By 
calculating the energy of such collective excitations (i.e. 
the matrix element (Flbq^H Landau bl, \F), where IF) is 
the electronic Fermi ground state defined in Appendix 
IB|) . it is straightforward to find in such subspace H Landau 

(CI) 

can be replaced by ^ ftwo blbq^ , where we have omit- 



ted a constant energy term. 

In order to express Hint in terms of second quantized 
operators, we have to calculate the matrix elements of 
the form 



{n,k,m\p-A,,n{r)\n',k',m') cx ^ {m\sm{nn,z/ L,)\m') {n\py\n') (a,,,„.(fc|e*«-"|fc') + aL,n,(fc|e"''=="|A:')) , (C4) 

I 



q^-n^ 



where |n, k,m) is the state of an electron in the m*'' hand side of Eq. (jC4p that m — m' . 
subband and n^^ Landau level with momentum k (see 
appendix [B|) . As already mentioned if we assume that 

Lqw < Lz, it comes from the first term in the right The matrix elements for the transition between Landau 

levels read 
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{n\py\n') = I™ ^° (Vn' + 1 J„,„'+i + \/n' (5„,„'_i) 
V2 (n — n') ^ ' 

Therefore, we get the harmonic oscihator selection rule 
n' = n ± 1; transitions such that |n' — n| > 2 are strictly 
forbidden. 

The integration over the x direction gives rise to the 
momentum conservation and we obtain the result 

qx,k 

^ ^g^"?^ \^''i^,/c-g^''!y-l,fc '^v-l,k'^v,k+q^J ■ 

The coupling constant (we are using the same notation 
as in section Hi D[) is given by 

I 



Given that the quantum well size Lqw ^ Lz, the 
integral along z gives the same selection rule as for 
Hint '■ m — m' so {m\siii('Klzz/ Lz) sin^nV^z/ Lz)\m') « 
sin(7rnz/2) sm{Trn'j2)6m,m' ■ 

On the other hand, the matrix elements along the other 
directions x and y can be factorized. For example the 
term (n, fc|e*(''^^'°')'^|n', fc') can be factorized as 

(fc|e<«^-«-)nfc')(n,fc|n',fc'). (C8) 

The second term of the previous equation is then given 
by the overlap integral (k, k') 

In,n' (k, k') = dy (j)n (y - yo(fc)) 0n' {v - Voik')) ■ 
Jo 

(C9) 

Assuming Ly » Iq (which is realistic being Iq sa 
O.l/Ltm for Wo = 52GHzrad~^) we get : 

/+CXD 
dy 4>n (y) <^«' iv - Aj/o) • (CIO) 
-oo 



Introducing the bosonic bright mode creation operator 
as described in section III Dl 



V nQWp2DEGS ^ 



we can check that the interaction Hamiltonian Hint 
can be exactly rewritten as it is given in Eq. (1101) by 
taking 



nQwP2DEGS 

V — I — 



Analogously in order to calculate Hdia we need the 
matrix elements : 



(C7) 



with Ayo = yo{k') — yo{k). We can then use the useful 
expansion of the harmonic oscillator wavefunctiong'^'^ 

0„ (y - Ayo) = E ('^' ) ^'^^ (y) . 

(Gil) 

where /3 = and L"^~" is the Laguerre polynomial 

of degree n and index m — n. Substituting Eq. (jGlip 
into Eq. (|C10|) . we obtain 

77M / 1 \ n—n , , 

/n,„' (fc, fc') « _ (/32) l;:,-" (/3) . (G12) 

Because of the gaussian dependence oc e "'o , 
-fn,n' (fc, k') exhibits a peak centered at Ayo = as sharp 
as /q- Thanks to the definition yo{k) = l^k and the 
condition that << Ly, we can consider approxima- 
tively that /„,„' {k,k') oc S^^k' ■ In addition, due to the 
orthonormality of the function 0„, setting k = k' (or 



{n, k,Tn\Al^\n',k',m') oc J2q^,q'^,n,,n'JM sin(7rn^z/L^) sin{Tm'^z/Lz)\m') 

+«L.i. 4,,,, (n,fc|e-<«^+«->|n',fc') a,, (n, fc|e<''^+«-)nn', fc') } • 

I 
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equivalently Ayo = 0) into the integral (IC10|) implies 

that /n,n' {k, fc ) ~ Sn ni. 

Hence, setting k = k' into the integral given by the first 
term of eq. (|C8[) . we obtain the corresponding selection 



rules on and q'^. Note that in the case of the previous 
example, we get that q'^ = qx- Finally, we find that H^ua 
is given by Eq. ([T0|) with the corresponding coupling 
constants. 
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